Irreversible magnetization switching using surface acoustic waves 
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An analytical and numerical approach is developped to pinpoint the optimal experimental con- 
ditions to irreversibly switch magnetization using surface acoustic waves (SAWs). The layers are 
magnetized perpendicular to the plane and two switching mechanisms are considered. In preces- 
sional switching, a small in-plane field initially tilts the magnetization and the passage of the SAW 
modifies the magnetic anisotropy parameters through inverse magneto-striction, which triggers pre- 
cession, and eventually reversal. Using the micromagnetic parameters of a fully characterized layer 
of the magnetic semiconductor (Ga,Mn)(As,P), we then show that there is a large window of acces- 
sible experimental conditions (SAW amplitude/wave-vector, field amplitude/orientation) allowing 
irreversible switching. As this is a resonant process, the influence of the detuning of the SAW fre- 
quency to the magnetic system's eigenfrequency is also explored. Finally, another - non-resonant - 
switching mechanism is briefly contemplated, and found to be applicable to (Ga,Mn)(As,P): SAW- 
assisted domain nucleation. In this case, a small perpendicular field is applied opposite the initial 
magnetization and the passage of the SAW lowers the domain nucleation barrier. 

PACS numbers: 73.50.Rb, 75.60.Jk,75.78.-n,75.50.Pp,62.65.+ik 



I. INTRODUCTION 

In a large number of ferromagnets, the coupling be- 
tween strain and magnetization originates from the spin- 
orbit interaction, and was shown early on to be maxi- 
mum when elastic and magnetic resonance (precession) 
frequencies match 1 -. This effect has been revisited in the 
light of spintronics applications in the past few years 
with compelling dynamic experiments in both magnetic 
semiconductors^ and metals 4 -. A first approach relies on 
the generation of picosecond acoustic pulses (longitudi- 
nal or transverse phonons). When coupled to the layer's 
magnons, magnetization precession may be triggered 3 , 
but it remains a fairly inefficient mechanism as the strain 
spectrum peaks quite high (20-30 GHz£) above typi- 
cal precession frequencies (0.5-10 GHz). Switching of 
a perpendicularly magnetized (Ga,Mn)(As,P) structure 
has recently been demonstrated using this technique 6 , 
but the effect was shown to originate from incoherent 
phonons (heat waves), and not from a magneto-strictive 
effect due to the high frequency coherent phonons pro- 
duced. Another route consists in generating strain 
through lower frequency (<2 GHz) surface acoustic waves 
(SAWs). On in-plane magnetized systems, SAWs have 
been used to drive ferromagnetic resonance in thin Ni 
films^, or periodically switch magnetization between hard 
and easy axes in Co bars 7 -. Recent theoretical work has 
focused on the switching of in-plane Terfenol nanomag- 
nets subjected to stress^, but no experimental or theo- 



retical work has been shown on perpendicularly magne- 
tized systems. These materials are for instance particu- 
larly relevant to high density magnetic information stor- 
age technologies. We believe SAWs offer two main advan- 
tages for magnetization reversal compared to picosecond 
acoustics: their relatively low frequencies, easily matched 
to precession frequencies, and the narrow bandwidth of 
the generated acoustic wave (a few MHz), as opposed to 
the broad-band spectrum in the former technique. 



In this work, we wish to address theoretically the irre- 
versible magnetization reversal in perpendicularly mag- 
netized layers using surface acoustic waves, and under 
realistic experimental conditions on a test system consist- 
ing in thin (Ga,Mn)(As,P) layers, a magneto-strictive di- 
lute magnetic semiconductor. Two possible mechanisms 
are considered, both relying on the transient modification 
of the magnetic anisotropy by the SAW. In precessional 
switching, the magnetization is pulled away from equilib- 
rium by an in-plane field, and the SAW triggers a large 
angle precession of the magnetization which may end up 
in a full reversal. In SAW-assisted domain nucleation, 
a small perpendicular field is applied opposite the ini- 
tial magnetization, and the SAW is used to locally lower 
the domain wall (DW) energy, and thus initiate domain 
nucleation, leading to a full reversal. 
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FIG. 1. (a) Set-up geometry for a SAW propagating along 
[100], and coordinates, (b) Depth dependence of the ampli- 
tude of the Rayleigh wave components (f—1 GHz) plotted 
using Annex B equations. 



II. DESCRIPTION OF THE SYSTEM 

1. Generation of SAWs 

SAWs are excited and detected by interdigital trans- 
ducers (IDTs) on a piezoelectric laye r 10 ' 11 deposited on a 
magnetic thin film (Fig. QJ,). We will for now limit our- 
selves to the case of a Rayleigh wave propagating along 
the [100] axis of a cubic crystal. The case of a wave prop- 
agating along [110] will be discussed in Section IV. 2. The 
only finite propagating strain wave components are then 
e xx {x, z, t), e zz {x, z,t), and e xz (x,z,t) (details in Annex 
B, axes defined in Fig. [IJi). Their wavelength is given di- 
rectly by the IDT period A R « 3-5 /im for /=0.5-l GHz, 
and their dispersion-free velocity by the elastic constants 
of the material, Vr=2711 m.s~\ with A R =^F. Their 
depth-dependence is plotted in Fig. [Tp. Two hypotheses 
may then be made if the magnetic layer is taken much 
thinner (<50 nm) than A^: (i) the e xz component can be 
neglected, as its amplitude remains weak close to the sur- 
face, and (ii) the strain field can be considered constant 
along z. We will therefore take z=0 in the expressions 
of e zz , e xx and e xz —0. Finally, the RF power passing 
through the combs is small enough (10 mW) to neglect 
any resulting heating of the sample. 



2. Magnetic system 

The time-dependent dynamics of the magnetization 
are described by the Landau-Lifshitz-Gilbert (LLG) 
equation: 



dM 
~dt 

VoHeff 



-jM X voHeff 
-VmF(M) 



a „-. dM 

M- Mx ln 



(i) 

(2) 



where M[9(x, t), 0(x, t)} is the magnetization expressed 
in polar coordinates with M s its norm (taken to be 

constant) and 7 > the gyromagnetic ratio. H e ff is 
the effective field, i.e the sum of a magneto-crystalline 
anisotropy term, a shape anisotropy term, and finally the 
Zeeman contribution from the externally applied field. 
In this work (except in Section V), the exchange energy 
contribution will be neglected and we will effectively be 
looking at the behavior of a single macrospin. 

Following Linnik et. al^, a normalized free energy 
density Fm—F/M s is defined, where F is a very general 
form of the free energy density of a cubic ferromagnetic 
layer distorted by strain: 

F M (M) = (A 2e - 2A 4s )e(x, t) cos 2 9+ 
(B c + 2A 4£ e(a;,t))cos 4 6+ 

isin 4 6»(B c - Ai e £(a;,t))(3 + cos40)+ (3) 

UqM s o 1 2 

— - — cos + -A 2xv e X y sin 9 sin 2<j>- 

[sin 9(hqH x cos <f> + fioH y sin 4>) + ^oH z cos 9] 

H e xt={H x ,Hy,H z ) is the externally applied field and 
9{x,t) 1 (j){x,t) is abbreviated into 9,(f>- B c is the cu- 
bic anisotropy field, and A^, A 2e ^ A 2xy are the mag- 
netoelastic coefficients. The dependence of the mag- 
netic anisotropy on strain is given through the terms s xy 
(static shear strain, details in Annex A), and s(x,t)— 
Aeo + Se(x, t). Ss(x, t)—Se zz (x, t) — 5e xx (x 1 1) is the strain 
generated by the SAW. The Rayleigh wave propagating 
in a cubic material along <f//[100] has been calculated 
analytically, and found to be quite different from the one 
used for isotropic materials [Ref JYQ] (details in Annex 
B). The resultant strain difference is given by 13 : 



5e(x, t) 



■, cos(o;i — qx) 



(4) 



where e max is the SAW amplitude 1 ^, and q its wave- 
vector, with q=uj/VR. Aeo=e zz ,o — e XX fi is the differ- 
ence between the static out-of-plane and in-plane strain 
components, related by the elastic modules Cn, C12 of 
GaAs 1 ^ through e zz , =-2§^e xxfi . 



III. ANALYTICAL SOLUTION: SMALL ANGLE 
PRECESSION 

1. Principles of precessional switching 

In so-called precessional switching, the perpendicular 
magnetization M is first tilted towards the layer by 
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an in-plane magnetic field. A short perturbation (e.g 
an optica l 15 ' 16 , acoustic 3 , or ultra-fast magnetic^! or 
electrioAS, field pulse) then modifies the micromagnetic 
parameters enough to change the effective field seen by 
the magnetization, and send it precessing. If the pre- 
cession amplitude is sufficiently large, the magnetization 
can switch to another potential valley, where it will re- 
main if the perturbation lasts an odd multiple of half the 
precession period^, or if damping eventually prevents 
M from oscillating between the two minima ("ringing" 
phenomenon). This mechanism has for instance being 
suggested for micro-wave assisted switching at a head 
field significantly below the medium coercivityi^ or for 
subnanosecond spin torque switching in magnetic tunnel 
junctions 20 . 



The arrival of the SAW on the magnetic layer leads to 
a modification of the magneto-strictive anisotropy terms, 
and thereby of the effective field H e ff. The time re- 
sponse of the magnetization is assumed short (around 
100 ps) on the SAW's time-scale, as evidenced by recent 
pump-probe experiments 2 .. Eq. ((T|) then shows that this 
triggers magnetization precession as long as the torque 
M x /io-ffe/ / remains non-zero and the damping has not 
aligned M back along the applied field. Two computa- 
tional approaches were then followed. Firstly, the con- 
ditions leading to magnetization precession were estab- 
lished by assuming small changes in magnetization direc- 
tion, 59, 6<p, in order to solve this equation analytically. 
Secondly, in view of establishing the experimental con- 
ditions leading to irreversible precessional switching of a 
(Ga,Mn)(As,P) layer, the LLG equation was solved nu- 
merically, and a switching diagram established. In this 
work, perpendicularly magnetized layers were considered. 
This is often a problematic configuration, since the en- 
ergy barriers are high for a full 7r reversal of the magne- 
tization. 

I 



2. General solution 

In this first approach, perturbations are small, leading 
to small changes in the magnetization direction around 
its equilibrium position Mo[0o,(f>o\. Provided the mag- 
netic anisotropy and applied fields are such that 9 ^0, 
Eq. ([T]) can be linearized into: 



-56 = 
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[F^S4> + F^gSB + F^Se] + ad'<p sin O (5) 
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- FeeSe] 
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50 (6) 



The terms Ft. 



stand for |^-. 



and the dot denotes 



the time derivative. In the following, the magnetization 
precession amplitude 50{x,t) will be calculated in x=0, 
but can easily be obtained at any distance x from the 
comb by computing 50(t — x/Vr). 

The eigen-frequency of the system in the absence of 
acoustic wave is first determined by assuming harmonic 
solutions for the angle deviations: 50 = 50$e 



f> Q e 



;n P t 



Expressing the determinant of the correspond- 
ing coupled equations system (5,6) then yields the com- 
plex precession frequency f2p of the magnetization in 
the presence of a finite damping term where we define 
Qp=ujp+ix'- 



2 2JJ2 



LUp = 



w 



VI + a 2 ' 

2(1 + a 2 ) 

7 

sin0o 
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F 2 



(7) 

(8) 
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We have further defined an effective field H a =Fgg + 
F§§l sin 2 0o, and the precession rate luq in the absence of 
damping. Assuming the SAW arrives at an instant t=0, 
the variation-of-parameters method then yields the am- 
plitude of the magnetization precession 50 (t) as a func- 
tion of the exciting SAW frequency u> and amplitude 
Smaxi the precession frequency wp and the damping a: 



50 it) = 



E-max^B 



(i 



J 2 + r 4 



f(u,/3) cos(ut + r]) - 



u; e 



-xt 



: cos(o;pf + £) 



(10) 



where £ and r\ are two phase shifts that depend on cj 
and the material's parameters. /(w,/3), j3 and fig are 
defined in Annex D. We further define: 



2 2 2 

U res = W P — X 



\/ 2uJ PX 



(11) 
(12) 



ui res is the resonance frequency of the system and T 
is related to the resonance broadening. This very gen- 
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eral expression of the precession amplitude highlights two 
physical behaviors. The first one is that, as expected in- 
tuitively, the precession consists of a forced term oscil- 
lating at the excitation frequency w, and a damped term 
at the eigenfrequency of the magnetic system uip. The 
second one if that the excitation frequency giving the 
largest amplitude is not exactly wp, but a slightly lower 
value, uj res , which is a modified resonance frequency of 
the damped system in the small perturbation regime. Fi- 
nally, a broadening term T allows the precession ampli- 
tude not to diverge at resonance. 

The amplitude of the precession is linear in e max and in 
fig which depends non-trivially on the magneto-strictivc 
coefficients, the damping, and the applied field through 
the value of 9q, 4>o (details in Annex D). We will see below 
that its expression can however be greatly simplified in 
some limiting cases. 



3. Application to thin (Ga,Mn)(As,P) layers 

The dilute magnetic semiconductor (Ga,Mn)(As,P) is 
a good test-bench material to investigate fast acoustics- 
induced magnetization switching. The carrier-mediated 
nature of its ferromagnetic phase results in a strong de- 
pendence of the magnetic anisotropy on the strain state 
of the layer, through the band splitting it induces on 
the light- and heavy-hole bands^i. For instance, Glunk 
et. al 2 ^ have shown that the perpendicular uniaxial 
anisotropy term in (Ga,Mn)(As,P) is proportional to 
both the out-of-plane strain coefficient e zz and the hole 
concentration p. Moreover, contrary to metals, typical 
precession frequencies of (Ga,Mn)(As,P) can be fairly 
low, of the order of the GHz in small magnetic fields 5 , 
which allows good matching to the acoustic wave frequen- 
cies provided by SAWs. Finally, the damping parameter 
can be rather high in this material (a=0.1-0.3 23 ' 24 ) com- 
pared to metals (a=0.01 in Ni 80 Fe2o), which will limit 
ringing effects preventing irreversible switching. Whereas 
in metals, precessional switching is mainly governed by 
the precession of the magnetization around the demagne- 
tizing field, in (Ga,Mn)(As,P) this process will be driven 
by the magneto-crystalline anisotropy since its magneti- 
zation at saturation is weak. 

While Eq. ([3]) conveniently highlights the magneto- 
strictive terms, the following form of energy is more com- 
monly used^ 5 - to determine experimentally the anisotropy 
coefficients in (Ga,Mn)(As,P) (details in Annex A for the 



correspondence between both energy forms): 

F M (9, 0) = -B 2± cos 2 9 - ^B 4± cos 4 9- 

1 7T 

-B if/ sin 4 6»(3 + cos 40) - B 2 1 1 sin 2 9 sin 2 (0 — —)+ 

— — - cos 2 9 — [sin9(iioH x cos + Mo-ffy sin 0) + /ioH z cos 9] 

(13) 

The magnetic anisotropy is largely dominated by the 
uniaxial term B 2 ±, followed by the cubic terms B 4 ± and 
B^j j which result from the tetragonal distortion of the 
lattice as the magnetic layer grows strained upon its sub- 
strate. A linear dependence of the uniaxial anisotropy on 
strain has indeed been found experimentally using vari- 
ous technique o 22 ' 26 ' 27 . The in-plane uniaxial term B 2 // is 
weakest and corresponds to a minor anisotropy between 
[110] and [110] axes (details in Annex A). 

To estimate quantitatively the amplitude of the pre- 
cession, a sample of relatively small perpendicular 
anisotropy is chosen in order to have GHz or sub- 
GHz precession frequencies, adapted to SAWs excited 
by micron-wide IDTs. An existing 50 nm thick sam- 
ple, with a static strain e 22i o=-0.05% and XMn ~ 7% 
is considered. This small lattice mismatch, yielding a 
moderate magnetic anistropy, is obtained by co-doping 
the (Ga,Mn) (As, P) layer with Phosphorus (yp ~ 4%) as 
described in Ref. H|| At 95 K, M s = 9 kAm" 1 and ferro- 
magnetic resonance spectroscopy yields: B 2 ± = 22.5 mT, 
B 4± = -2.3 mT, B 4// = 2.3 mT and B 2// = -1.2 mT. 
The damping will be taken as a=0.1, but note that this 
term has been shown to vary between 0.001 and 0.3 with 
magnetic and electric doping, as well as whether one 
measures the extrinsic damping or an intrinsic Gilbert 
dampin g 23 ' 24 ' 29 . 

Let us first put some numbers on the relevant fre- 
quencies (/fe=Wfe/27r). Under an in-plane magnetic field 
of 2 mT (6»o=3°), Eqs. (7,9) yield: / =1.017 GHz, 
./p=1.007 GHz, / res =1.002 GHz. The decrease of the 
resonance frequency due to the inclusion of damping is 
therefore relatively small, a mere 1.5%. The broadening 
is rather average: 90 MHz (full-width at half maximum). 
Finally, the exponential damping of the precession occurs 
on a time-scale of l/x=1.6 ns. 

In order to isolate the relevant parameters to obtain a 
large angle precession, ttg and /(o>,/3) may be simplified 
provided the explicit energy density of (Ga,Mn)(As,P) 
(Eq. (fTB")) ) is used and a few hypotheses are made. Since 
in general /3 < 1 (see Annex D), we develop /(w,/3) w 
w p + X0- We also use: A ie « A 2 J£, a « 1, and 
consider that 0o closely follows the applied field direction 
4> H , so that fi H x cos O + HoH y sin O « ^H ext . The 
precession amplitude at resonance can then be simplified 
into: 



\se\, 



' p2 



9,0 



8f 
w 



A 2e cos 9 (B 4 // sin 3 6> cos 40 o — B 2 // sin 6> sin 20 o + ^oH ext /2) (14) 



In Eq. (|T4l) . ^4 2 e is roughly proportionnal to the uni- 
axial anisotropy term B 2 ± (see Annex A) while the in- 
plane anisotropy terms B^i / and B 2 // are affine func- 
tions of A 2e . One can see that a larger precession am- 
plitude | (5$ | mQ;r first requires a large uniaxial anisotropy, 
and large in-plane anisotropics -B4// and B 2 //] but this 
will tend to increase precession frequencies high above 
typical SAW frequencies for micron-sized IDTs. The ap- 
plied field amplitude and angle can however also be opti- 
mized as we show in the following numerical calculations 
of |^0|max at fixed strain amplitude e ma x=10 -5 (Fig. [5]). 
This value is taken deliberately small to remain in the 
small perturbation regime. 

At fixed field amplitude /ioH ext =24 mT (large enough 
to insure 0o ~ 4>Hi $o=45°), the angle of the field is 
first varied in the plane (Fig. [2t0- The precession ampli- 
tude is largest in the (/>#=0-90° range, with a maximum 
at </>#=45°. This results from the competition between 
the two in-plane anisotropics terms maximized at 0o=O° 
modulo 90° (for S 4// > 0), or at </> =45° (for B 2// < 0). 
The amplitude variations are however weak, and this is 
clearly not the most critical parameter. At fixed field 
angle (f>H=0° this time, the amplitude of the field (and 
therefore of the initial tilt 9 ) is made to vary (Fig. Od). 
The variations observed are this time more pronounced, 
and 1 59 \ max is clearly maximum when the magnetization 
is most pulled away from its zero-field orientation, as was 
also concluded from theoretical studies of precession trig- 
gered by pico-second acoustic pulses 1 ^. For fields above 
30 mT, the magnetization is saturated in the plane, and 
the precession amplitude plummets down to a few 10 -4 
rad, and gradually decreases to zero. The precise varia- 
tion of precession amplitude with field amplitude and ori- 
entation of course depends on the value of the anisotropy 
parameters, but since up, /3, and V vary slowly with these 
parameters, the analytical dependence of fls,o with (fin 
and fioH ext given in Eq. (fT4l) gives a good idea of the 
conditions maximizing this amplitude. 

Finally, note that we have not taken into account 
the influence of the ferromagnetic resonance upon the 
acoustic wave propagation. We have indeed assumed the 
phonon-magnon coupling in (Ga,Mn)(As,P) sufficiently 
weak to neglect in first approximation the absorption of 
the acoustic wave upon interaction with the ferromag- 
netic layer. Please refer to Ref. [3(3 for a complete ana- 
lytical treatment of this so-called "back-action" effect. 

In summary, analytically solving the LLG equation in 
the presence of a SAW of given frequency has allowed 
us to identify the experimental parameters apt to yield 
the largest precession amplitude in a perturbative regime: 
field as large as possible without saturating the layer, and 
applied between [100] and [010] axes. 




(a) M de 9) ^o H ex,( mT ) (b) 

FIG. 2. Conditions maximizing the precession amplitude 
in a perpendicularly magnetized (Ga,Mn)(As,P) layer for 
£max=W~ 5 . (a) At fixed field amplitude, variation of the 
precession amplitude as a function of the field angle <f)H- (b) 
At fixed field angle <j>H~0° , variation of the precession ampli- 
tude as a function of the field amplitude /j,oH ext (red symbols) 
compared to the tilt 80 (solid black line) before the arrival of 
the SAW. 



IV. NUMERICAL SOLUTION: IRREVERSIBLE 
PRECESSIONAL SWITCHING 

1. Conditions for precessional switching 

To explore the conditions for precessional switching, it 
is necessary to go beyond the small angle approximation 
and solve Eq. (fTJ) numerically. The same sample as above 
is considered, with the field applied in the plane of the 
layer, along (/>#=0°. 

Experimentally, SAWs can be excited by rf bursts of 
length Tsaw ~ 150ns 10 . Here, we will only consider 
what happens during a single period. The rise-time r 
of the signal is given directly by the transit time of the 
acoustic wave through the emitting IDT. For about 10 
pairs of teeth working at sub-GHz frequencies, a realis- 
tic value is r w 20 ns. The experimental time profile 
of e(x, t)=A£o+Se(x, t) taking into account a rise and 
decay time (linear experimentally, but modeled as expo- 
nential) is computed as shown in Fig. [3^. The SAWs 
line-width is roughly given by 1/Tsaw— 7 MHz. Four 
main parameters can then be adjusted numerically to ex- 
plore the different behaviors of the system: the SAW am- 
plitude, e max (5.10 _5 -10~ 3 ), the in-plane magnetic field 
amplitude, which in turn controls the initial tilt of the 
magnetization 00, the detuning of the SAW frequency to 
the precession frequency , and finally, the damp- 

ing parameter a. Of those parameters, the first two can 
easily be changed during an experiment. Note that the 
detuning can equally be defined with respect to f re s, as 
they are within 1% of each other. 

At fixed SAW amplitude and damping, two distinct 
behaviors are observed. Examples are shown in Figs. 
|3}j,c, where the magnetization initially points upwards 
(M z /M s =l) before the application of the field, and the 
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(a) time (ns) (b) time (ns) (c) time (ns) 

FIG. 3. Time behavior of the magnetization from numerical simulations for £ maa; =2.5.10 -4 , a=0.1, Sf/fp=5% and an initial 
magnetization pointing towards the upper half, M z /M s ~ 1. The field is applied along [100]. (a) Temporal profile of e(x=0,t) 
for an rf burst Tsaw = 150ns at is aw = 700 MHz, and a rise time r=20ns. (b) Large angle precession, (c) Sustained switching 
leading to an irreversible reversal at the end of the SAW excitation. The SAW excitation time is indicated by the red dashed 
line. 



precession frequency lies around fp ss 780-930 MHz for 
the fields investigated. At low field (11 mT), such that 
the initial magnetization is moderately tilted towards 
the plane (9o=18°, Fig. the magnetization remains 
pointing up during the excitation, and precesses at is AW 
in a cone that is widest when the SAW has reached its 
stationary regime. At the SAW extinction, the magneti- 
zation returns to its initial position. This regime will 
defined as "large angle precession". Indeed, the am- 
plitude of this precession is about 10-100 times larger 
than the one observed in picosecond acoustics-triggered 
precession^: the strain pulse amplitude at the precession 
frequency is weak in this latter technique. At larger ap- 
plied field (18 mT), the magnetization first precesses in 
the upper quadrant at isAW, fully switches to M z /M 8 s» - 
1, and then oscillates between up and down positions 
at half the excitation frequency (Fig. By adjust- 

ing T S aw in the range [2n/fsAW,2{ n + tf/fsAw], the 
magnetization can then be released in the downwards 
position. This regime will be named "sustained switch- 
ing, conditional reversal" . It relies on the same mecha- 
nism as the precessional switching using tailored ultra- 
fast magnetic field pulse which was used in garnet o 17 ! 31 
but the lower precession frequencies make the adjustment 
of Tsaw less constraining. The optimal value of Tsaw 
for reversal is quite experiment-dependent though, as it 
depends on both the decay time of the SAW (related to 
the number of teeth in the IDT), and to the magnetiza- 
tion damping a. 

In between those two behaviors is a transition regime 
(not shown) where the magnetization undergoes quite a 
chaotic behavior, at times switching irreversibly before 
the end of the SAW. Finally, for an applied field saturat- 
ing the magnetization in the plane of the sample (9o—90°, 
not shown), the SAW induces a precession of M around 
the applied field, with an amplitude that decreases when 
the strength of the applied field increases. 

A more thorough exploration of the parameter space 
is shown in Fig. Hh,- The SAW amplitude e max is varied 
in steps of 2.10~ 4 or 2.10 -5 , the in-plane applied field in 



steps of 1 mT with </>#=0°, and the frequency detuning 
is first fixed to 5%. The resulting diagram shows that the 
behavior of Fig. [3] is very generic, regardless of the SAW 
amplitude: a large amplitude precession regime at low 
fields (in black), a conditionnal switching regime at high 
fields (in gray), and in between a transition regime (white 
line). As expected intuitively when e max decreases, the 
field necessary to obtain switching increases in order 
to keep the precession wide and compensate for this 
lesser efficiency Note that for large strain amplitudes 
{smax >5.10~ 4 ), this generic behavior is maintained, but 
multiple frequencies appear in the sustained switching 
regime due to strong non-linearities. The critical fields 
obtained for a larger detuning (20%) are indicated by the 
dashed line. Quite counter-intuitively, these seem to be 
lower than the ones found for <5f/f=5%, in particular at 
higher strain amplitudes, as if a larger precession ampli- 
tude were obtained away from resonance as opposed to 
at resonance. To elucidate this, we set the field to 11 mT 
and systematically recorded the precession amplitude as 
a function of the frequency detuning, for three differ- 
ent strain amplitudes. The result is shown in Fig. |4Jd, 
and confirms that whereas at low strain amplitude, the 
maximum precession amplitude is indeed obtained at res- 
onance (<5f/f=0%), when the strain amplitude increases, 
the maximum precession amplitude can be obtained quite 
far from resonance, at an increasingly large frequency de- 
tuning. This confirms that the small-perturbations ap- 
proach of Section III is only valid for small strain am- 
plitudes, and that beyond, the behavior becomes highly 
non-linear and the system's eigenfrequency is very prob- 
ably not given by fp anymore, but by a lower frequency. 
Finally, the main effect of decreasing the damping (not 
shown) is to lower the critical field between large angle 
precession and precessional switching. The important 
conclusion of these simulations is that there is a large 
region of the {e max , ^oH ext ) parameter space where ir- 
reversible switching of a macrospin is possible, and this 
at fairly low fields. 

Note that a more elegant way of identifying switch- 
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ing conditions would be to determine an analytical cri- 
terion leading to the growth of non-linearities in the sys- 
tem. This is not trivial in this coupled 9=f(<p,9,9), 
4>=g{9, <fi, 4>) system, and is beyond the scope of this pa- 
per. 



2. Particular cases: i) SAWs propagating along (110), ii) 
buried layers 

Here we address the case of a Rayleigh wave prop- 
agating along a (110) direction, a configuration easier 
to implement experimentally. A straightforward calcu- 
lation using a 7r/4 rotated frame shows that the SAW 
once again has three components, denoted Exx, £zz 
and exz, of identical shape and similar amplitude as 
the one traveling along [100] (details in Annex B). Given 
this strain tensor, we rotate it back into the x//[100] 
frame and inject the resulting components in the energy. 
The free energy density of the layer is then identical to 
Eq. (O with e(x, t)= Ae + e zz {x, t) — e xx (x, t) replaced 
by: Ae + £zz(x,t) ~ ^s X x{x,t), and e xyfl replaced by 
e xy , ± \e X x (+ for g//[110], - for <?//[110]). The cal- 
culation therefore requires knowing e XVt o which is prob- 
lematic: there is no experimental measurement of this 
shear strain, and therefore no independent determination 
of A2 xy and e xy .o is possible. We will therefore stick to 
the £a;y,o=10~ 4 value used in Ref. [H, and deduce from 
our experimental value B 2 // = -1.2 mT the parameter 
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The amplitude of the precession triggered by the SAW 
is then estimated numerically using the same micro- 
magnetic parameters as above. For strain amplitudes 
£maa;=0.5 to 10. 10~ 4 , and /ioH ext =5 mT along [100], the 
precession amplitude when the SAW is traveling along 
[110] or [110] is systematically about half of when it trav- 
els along [100]. This is mainly due to the fact that the 
strain amplitude in front of the A% s coefficient is reduced: 
the biaxial strain is along the crystallographic axes (100), 
and effectively reduced along the diagonals (110). The 
magneto-strictive process is then less efficient. 



Finally, we look at what happens for either a thick 
magnetic layer, or a buried layer, since the three Rayleigh 
wave strain components oscillate in amplitude with z 
(Fig. Qi>). Away from the surface, e zz and e xx both re- 
duce to their first zero within about 0.1-0.2Asaiv, whilst 
e xz takes over, albeit with a much smaller amplitude. For 
layers of this order of thickness (<300 nm in our case), 
the conclusions drawn above remain valid. However, an 
interesting case is that of a thin magnetic layer buried 
where e xz is maximum, and 6s non-zero. As evidenced 
in Fig. [lb, this lies around z c =0.25Asaw, where Se is 
then still at about half its surface value. 

In the (Ga,Mn)(As,P) free energy density form (Eq 
(|13p ). there is no dependence on e xz as this strain com- 
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FIG. 4. (a) Switching diagram of a macrospin excited by a 
SAW with q=0.1, a SAW frequency detuned by 5% from the 
precession frequency fp, and an in-plane field along [100]. 
The large angle precession regime is separated from the sus- 
tained switching regime by the continuous white line (transi- 
tion regime). Critical fields for a frequency detuning of 20% 
are indicated by the dashed line, (b) Amplitude of the pre- 
cession £^(Fig. [3b behavior) as a function of the detuning, 
for different strain amplitudes. 



ponent is not present statically in the layers. However, 
just as a e xy term will give a magnetic anisotropy be- 
tween [110] and [110] axes (term \Ai xy e xy sin 2 #sin2</> 
in Eq. ([3])), it is easy to imagine that a e xz term will 
give a magnetic anisotropy between [110] and [001] axes, 
expressed as: \A2 XZ s xz sin2(9cos(/>. This term was eval- 
uated numerically to Ai xz =%§ T (details in Annex C) 
close to 95 K. 

The amplitude of the magnetization precession trig- 
gered by the SAW was then compared at fixed applied 
field 5 mT, and SAW frequency 2 GHz (5% detuned from 
the precession frequencies that are larger than at 95 K 
since the anisotropy is stronger) and strain amplitude 
£max—^0~ 5 for a surface layer (z=0) or for a layer buried 
at z c —0.25Asaw- In the latter case, the precession am- 
plitude is about three times larger compared to a surface 
layer. This is due to the fact that A2 XZ is large and that 
the three strain components now contribute to the pre- 
cession. However, for the precession frequencies explored 
at 95 K, As aw is too large (several microns), and z c 
therefore unrealistically deep. This option might how- 
ever prove more relevant for materials requiring higher 
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SAW frequencies (smaller Asaw)- 



V. SAW ASSISTED DOMAIN NUCLEATION 

Let us mention briefly another approach for irreversible 
switching induced by a SAW. More likely than the (co- 
herent) precessional switching of a structure, is the pos- 
sibility that the SAW will locally lower a domain nu- 
cleation barrier and thus switch the whole layer un- 
der a small (propagating) magnetic field applied con- 
currently opposite to the initial magnetization. Indeed, 
in perpendicularly magnetized layers, the coercive field 
HqHc is largely determined by DW nucleation and/or 
propagation barriers 32 . Transient reduction of coerciv- 
ity has already been demonstrated in garnets 3 - 3 , or mag- 
netic semiconductors 3 -^ 3 -^ using ultra-fast light pulses. It 
is also at the basis of thermo-magnetic writing. 

Neglecting the stray field energy between the nucleated 
domain and the rest of the layer—, the domain nucleation 
barrier can be expressed as E nuc =2Trr nuc d x er, where d 
is the layer thickness, r nuc the radius of the nucleated 
domain, and a its surface energy. The latter depends on 
the perpendicular uniaxial anisotropy and the exchange 
constant A ex through a cx yj 'A ex B± (e xx , e zz )3&: If the 
magnetization reversal is nucleation limited, switching 
then occurs on a typical timescale given by the Arrhe- 
nius law t—t exp( Ez ^^p' uc ) > where fc^T is the thermal 
energy, and Ez=-HoH ext M s x dirr 2 is the Zeeman energy 
which lowers this barrier, with H ext the field applied per- 
pendicular to the layer, tq is the typical time needed for 
an energy exchange between spin and lattice, and de- 
pends on the damping, the anisotropy constants and the 
magnetization 3 -^. It is usually estimated to about 10 ps. 

For the sample at 95 K considered above, the effective 
uniaxial anisotropy field is B_L«i?2±=23 mT. Previous 
temperature dependent experiments on (Ga,Mn,As,P) 38 
have moreover given A ex =10~ 13 pJ/m, providing an 
estimation of the DW surface energy cr=2.10~ 5 J/m 2 . 
Experimental hysteresis loops done at 90 K yield 
fj,oH c —2 mT, with the perpendicular field applied in 
«100 ms pulses. Using the Arrhenius law, this gives an 
experimental estimation of the nucleated domain's diam- 
eter of a few nanometers. The nucleation barrier is then 
of the order of 44 meV (much lower than the intrinsic 1 eV 
barrier estimated using DW nucleation theory- 2 ^), the 
Zeeman lowering around -3A0~ 2 H ext meV/mT, while 
the thermal energy lies around 8 meV. In the presence 
of the SAW, the effective anisotropy coefficient becomes: 
Bf AW {x,t)=(A 2e - 2^ 4e )(Aeo + Se(x,t)). Because the 
lattice mismatch is small in this sample (in order to 
have a weak anisotropy), the transient strain modifica- 
tion 5e(x,t) can very well be of the order of the static 
strain mismatch Ae (Fig. Ot), thus strongly reducing 
the uniaxial anisotropy field. This dramatically lowers 
the DW nucleation barrier during about a quarter of the 
SAWs period. If the resulting switching time r is shorter 
than the time during which this barrier is very low, ultra- 



fast magnetization switching is expected to occur under 
a small (propagating) perpendicular field. Note that this 
is a non-resonant process, and as such does not involve 
any strong constraint on the SAW frequency. On the 
contrary, it might be more relevant to aim for low fre- 
quency, which will leave plenty of time for the domain to 
nucleate during the transient decrease of anisotropy. 



VI. DISCUSSION 

In the first approach (Sections. Ill, IV), the calcu- 
lations and simulations were done for a macrospin at 
a fixed distance from the combs. For an extended, 
single-domain structure (no exchange energy), the con- 
clusions will remain identical as long as the shape 
anisotropy of the structure is not significantly different 
from the cos 2 9 term of Eqs. © , CE3]) . For micron- 
sized structures, and given the weak magnetization of 
(Ga,Mn)(As,P), this is legitimate. However, for samples 
with a strong shape anisotropy, a modified demagneti- 
zation factor would need to be included in the energy 
form, as was done by Roy et. al£. A more proper solu- 
tion would however eventually need to take into account 
the exchange contribution in the free energy. In both pi- 
cosecond acoustic pulses- and SAW-induced precession, 
the modification of the eigen-frequencies due to this ex- 
change contribution was shown to be negligible 3 ^, given 
the dispersion curve in the considered regime remains al- 
most flat. An estimation of this contribution to the real 
space, i.e the generation of spin waves to accommodate a 
spatially varying magnetization vector is a problem more 
appropriately tackled by micromagnetic finite elements 
methods, and is beyond the scope of this paper. 

In the second approach, a phenomenological model of 
SAW- assisted domain nucleation is used where exchange 
is implicitly taken into account (formation of a domain 
wall), and the field applied perpendicularly to the plane. 
Once again, SAWs seem more adequate than picosecond 
acoustic pulses. Indeed, were the transient strain strong 
enough at the relevant frequency to induce a substantial 
change of the uniaxial anisotropy coefficient to lower the 
domain's nucleation barrier, this would only last about 
10 ps^, which would require an ultra-fast DW nucleation. 
Looking at the Arrhenius law once again shows that this 
would mean effectively canceling the energy barrier, in 
order to have t«to, which seems unlikely given the weak 
accessible transient strain amplitudes in this technique. 



VII. CONCLUSION 

A general analytical approach to SAW-assisted magne- 
tization precessional switching has been developed based 
on the strain dependence of the magnetic anisotropy co- 
efficients, and taking the damping into account. Several 
parameters were found to be important, such as a large 
in-plane field, and SAWs propagating along [100] rather 
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than [110]. Numerical simulations using realistic exper- 
imental parameters of (Ga,Mn)(As,P) then clearly evi- 
denced a wide range of fields and SAW amplitudes under 
which irreversible switching was possible. SAWs were 
shown to possibly be a more adequate method to switch 
magnetization than picosecond acoustic pulses. Finally, 
although these concepts were tested on (Ga,Mn)(As,P) 
thin films, the analytical forms of energy and precession 
amplitude given in this work make them applicable to 
any magneto-strictive material. 

We gratefully acknowledge insightful advice from C. 
Tanguy at Orange labs, and technical helop from S. Ma- 
jrab and L. Becerra. This work was performed in the 
framework of the MANGAS project (ANR 2010-BLANC- 
0424-02). 

ANNEX 

A. Magnetic anisotropy coefficients 

The magnetic anisotropy terms A^, A 2e , A 2xy , and 
B c of Eq. (j3|), and the terms B 2 ±, B 4 ±, B 2 // and B4// 
(obtained experimentally by ferromagnetic resonance ex- 
periments for instance) of Eq. (TJ3]) are related as follows: 



s xx (r, z, t) = -2i£ q e 1 ^' 2 e~ a<lz cos(bqz + fa/2)e^ t ~^ 
e zz {r,z,t) — 2 i p ( a qe**°/ 2 e- a9 '[-awa.(bqz + fa/2 - fa) + 

bcos{bqz + fa/2 - fa)]e^ wt -^ 
e xz (r,z,t) = -Coqe l ^ /2 e- aqz [acos{bqz + fa/2) + 

bsm(bqz + fa/2) + psm(bqz + fa/2 - fa)]e i{ - ut ~^ 

We define the wave- vector q (norm q), the posi- 
tion r (norm r) and Co the amplitude of the dis- 
placement. When <f//(100), q.r=qx and for a mag- 
netic layer much thinner than Ksaw, we can set z=0 
(Fig. Q}, and simplify 5e(x,t)=e zz (x,0,t)-e xx (x,0,t) 
into: 8e(x,t)=e max cos(ujt — qx), where e max is the am- 
plitude of the resulting wave. 

The parameters are related by: a cos "02/2 + 
&sinV> 2 /2 = psin(fa - fa/2). For GaAs and <f//(100), 
the values found numerically are: fa=-0.328, fa=-1.9, 
,9=1.18, a=0.402 and 6=-0.561. For <f//(110), these co- 
efficients are only slightly modified: fa=-0.531, ^ 2 =-2.1, 
p=1.34, a=0.500 and 6=-0.480, yielding strain ampli- 
tudes about 15% smaller than along (100) directions. 



A 2 e 
B c 

A 2x y — 



B 4// - B 



4± 



6Ae 

B4// - B±± B 2 // — 2£? 2 j_ 



3Ae n 



2Ae n 



B x 



2B. 



B 



2// 



where Ae = e zzfi - e xxfi 



For the vast majority of (Ga,Mn)(As,P) samples, the 
relationship A 4e « A 2e holds. 

Note that there is no experimental evidence of an £ xy $ 
shear strain, but rather it is a physical effect of the 
same symmetry that is at the root of the weak uniax- 
ial anisotropy B 2 j /. During the growth, when atoms are 
mobile on the surface, nearest-neighbor Mn pairs on the 
GaAs (001) surface have a lower energy for the [1-10] 
direction compared to the [110]^. 



B. Derivation of the strain wave expression 



C. Effect of a e„ strain component on 
magnetization precession 



In the presence of a e xz strain component, one can 
expect a magneto-strictive component^ of the form 
A 2xz e xz m x m z —^A 2xz £ xz sin2#cos</>. This term can be 
evaluated theoretically using an effective mass Hamilto- 
nian with the six-band k.p Luttinger-Kohn term, a strain 
tensor, and the p-d exchange interaction of the holes and 
the Mn spins in the molecular-field approximatio n 21 ' 40 . 
The saturation magnetization was set to 6 kA.m -1 
(T sa 95 K), which corresponds to a \Bq\ ~4 meV spin 
splitting parameter—. The usual biaxial strain terms 
£zz,a,£ xx ,o were set to zero, and a non-zero term ids xz 
introduced in the Bir-Pikus strain tensor, where d—- 
4.8 eV is the shear deformation potential^ -. In this 
way, B c and A 2xz were the only unknown parameters 
in the free energy density. The energy difference F(9) — 
F([001]) = Bc(cos 4 (9 + sin 4 9) + \A 2xz e xz sin 29 was then 
computed and fit numerically for p=3.10 20 cm -3 yield- 
ing: ^2^^=80 T. This value is quite large, in fact larger 
than any of the anisotropy parameters, but the resulting 
anisotropy field expected to be less than 10 mT. 



The strain components of the Rayleigh wave propagat- 
ing in a cubic material along cf//[100] are quite different 
from the usual formulas found in textbooks for isotropic 
materials [Ref JYQ]. They can be found analytically as: 



D. Small angle magnetization precession amplitude 

The precession amplitude 89 given in Eq. (fTUf depends 
on various parameters given below: 
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fe 
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tan/3 



7 

sin 8 

7 



sin t) 



2 



2(1 + a 2 ) 



7 



2(l+a 2 ) 
sin #o 



(15) 
(16) 

(17) 

(18) 



[fya sin 6» - /op 



, sin 00(2^00 + a 2 {F u - F ee sin 2 0„)) + fe(aH a sin 2 <? + 2F e0 sinfl (l + a 2 ))] 5 



sin 4 o [(l + a 2 )^-a 2 # 2 ] 



(19) 



In the case of (Ga,Mn)(As,P), some approximations 
can be made using A± s < < A 2e and developing £lg to its 
zero-th order expansion in a: 



fe = 



-^A^ e sin 3 9o sm ' 
-27^4 2e cos Q 



n fl ,o = -i / e 2 + 



(jgjfy + /<A^ ) 2 7 2 
sin 2 # w 2 



(20) 
(21) 

(22) 
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